Quasinormal modes of Schwarzschild black holes in four and higher dimensions 



m 

O 

o 

(N 
D 

a^ 

> 



0^ ' 
O . 
m , 
o ■ 

I 



X 



Vitor CardosoQ Jose P. S. LemosQ and Shijun Yoshid^ 

Centra Multidisciplinar de Astrofisica - CENTRA, 
Departamento de Fisiea, Instituto Superior Tecnico, 
Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal, 
(Dated: February 4, 2008) 

We make a thorough investigation of the asymptotic quasinormal modes of the four and five- 
dimensional Schwarzschild black hole for scalar, electromagnetic and gravitational perturbations. 
Our numerical results give full support to all the analytical predictions by Motl and Neitzke, for 
the leading term. We also compute the first order corrections analytically, by extending to higher 
dimensions, previous work of Musiri and Siopsis, and find excellent agreement with the numerical 
results. For generic spacetime dimension number D the first-order corrections go as (p_3)''/(p_2) ■ 
This means that there is a more rapid convergence to the asymptotic value for the five dimensional 
case than for the four dimensional case, as we also show numerically. 

PACS numbers: 



I. INTRODUCTION 

The study of quasinormal modes of black holes began 
more than thirty years ago, when Vishveshwara no- 
ticed that the signal from a perturbed black hole is, for 
most of the time, an exponentially decaying ringing sig- 
nal. It turns out that the ringing frequency and damping 
timescale are characteristic of the black hole, depending 
only on its parameters (like the mass and angular mo- 
mentum). We call these characteristic oscillations the 
quasinormal modes (QNMs) and the associated frequen- 
cies are termed quasinormal frequencies (QN frequen- 
cies), because they are really not stationary perturba- 
tions. Not surprisingly, QNMs play an important role in 
the dynamics of black holes, and consequently in grav- 
itational wave physics. In fact, it is possible 0j 13 to 
extract the parameters of the black hole simply by ob- 
serving these QN frequencies, using for example gravita- 
tional wave detectors. The discovery that QNMs dom- 
inate the answer of a black hole to almost any exterior 
perturbation was followed by a great effort to find, nu- 
merically and analytically, the QN frequencies. For ex- 
cellent reviews on the status of QNMs, prior to 2000, we 
refer the reader to Kokkotas and Schmidt ^] and NoUert 

. It is important to note that on the astrophysical as- 
pect, the most important QN frequencies are the lowest 
ones, i.e., frequencies with smaller imaginary part, and 
the most important spacetimes are the asymptotically 
flat and perhaps now the asymptotically de Sitter. How- 
ever, three years ago Horowitz and Hubeny @ pointed 
out that QNMs of black holes in anti-de Sitter space have 
a different importance. According to the AdS/CFT cor- 
respondence |3 , a black hole in anti-de Sitter space may 



be viewed as a thermal state in the dual theory. Perturb- 
ing this black hole corresponds to perturbing the thermal 
state, and therefore the typical timescale of approach to 
thermal equilibrium (which is hard to compute directly 
in the dual theory) should as well be governed by the 
lowest QN frequencies. That this is indeed the case was 
proved by Birmigham, Sachs and Solodukhin Q for the 
BTZ black hole, taking advantage that this is one of the 
few spacetimes where one can compute exactly its QN 
frequencies, as showed by Cardoso and Lemos 0. A 
similar study was made by Kurita and Sakagami lu\ for 
the D-3 brane. This interpretation for the imaginary part 
of the QN frequencies in terms of timescales of approach 
to thermal equilibrium in the dual conformal field theory 
has motivated a generalized search for the quasinormal 
modes of different black holes in anti-de Sitter spacetime, 
over the last three years |^l9llll|. 

Recently, the motivation for studying QN modes of 
black holes has grown enormously with the conjectures 
111 III 111 111 relating the highly damped QNMs (i.e., 
QN frequencies with large imaginary parts) to black hole 
area quantization and to the Barbero Immirzi parameter 
appearing in Loop Quantum Gravity. The seeds of that 
idea were planted some time ago by Bekenstein A 
semi-classical reasoning of the conjecture that the 

black hole area spectrum is quantized leads to 



An ^ llpn , 1, 2, . 



(1) 



Here Ip is the Planck length and 7 is an undetermined 
constant. However, statistical physics arguments impose 
a constraint on 7 [isL |2(]| : 



7 = 4 log k , 



(2) 
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where k is an integer. The integer k was left undeter- 
mined (although there were some suggestions for it ^19|), 
until Hod 12], supported by some of Bekenstein's ideas, 
put forward the proposal to determine k via a version 
of Bohr's correspondence principle, in which one admits 
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that the real part of QN frequencies with a large imagi- 
nary part plays a fundamental role. It was seen numeri- 
cally by NoUert 0, |^ that QN frequencies with a large 
imaginary part behave, in the Schwarzschild geometry as 

wAf = 0.0437123 + ^ (271 + 1) , (3) 
8 

where M is the black hole mass, and n the mode number. 
Hod first reahzed that 0.0437123 and went on to 

say that, if one supposes that the emission of a quantum 
with frequency ^ corresponds to the least possible en- 
ergy a black hole can emit, then the change in surface 
area will be (using A = IBttA/^) 

A A = 3277 MdM = 32t: Mhuj = 4ft In 3 . (4) 

Comparing with we then get fc = 3 and therefore the 
area spectrum is fixed to 

An = 41n3/pn; n 1,2, ... (5) 

It was certainly a daring proposal to map 0.0437123 to 
In 3, and even more to use the QN frequencies to quan- 
tize the black hole area, by appealing to "Bohr's corre- 
spondence principle" . The risk paid off: recently, Dreyer 
[13] put forward the hypothesis that if one uses such a 
correspondence it is possible to fix a formerly free param- 
eter, the Barbero-Immirzi parameter, appearing in Loop 
Quantum Gravity. Dreyer's work implied that the en- 
tropy has contributions from J = 1 edges mainly. One 
interpretation (proposed by Dreyer), is that one changes 
the gauge group from SU{2), which is the most natural 
one, to 50(3). There arc other explanations, though. 
For example, Corichi ,15] presents an argument based on 
simple conservation principles to suggest that one can 
keep SU (2) and still have consistency with QNM, while 
Ling and Zhang 16] propose to consider the supersym- 
metric extension of the theory. 

All of these proposals and conjectures could be useless 
if one could not prove that the real part of the QN fre- 
quencies do approach i.e., that the number 0.0437123 
in NoUert 's paper was exactly This was accom- 

plished by Motl 22] some months ago, using an ingenious 
technique, working with the continued fraction represen- 
tation of the wave equation. Subsequently, Motl and 
Neitzke ^ used a more flexible and powerful approach, 
called the monodromy method, and were not only able 
to rederive the four dimensional Schwarzschild value , 
but also to compute the asymptotic value for generic D- 
dimensional Schwarzschild black holes thereby confirm- 
ing a previous conjecture by Kunstatter 24]. They have 
also predicted the form of the asymptotic QN frequencies 
for the Reissner-Nordstrom geometry, a prediction which 
was verified by Berti and Kokkotas (25»] . 

It is interesting to speculate why, in thirty years inves- 
tigating QNMs and QN frequencies, no one has ever been 
able to deduce analytically the asymptotic value of the 
QN frequencies, nor even for the Schwarzschild geometry, 
which is the simplest case. In our opinion, there are three 



important facts that may help explain this: first, the con- 
cern was mostly with the lowest QN frequencies, the ones 
with smaller imaginary part, since they are the most im- 
portant in the astrophysical context where QNMs were 
inserted until three years ago. In fact, most of the tech- 
niques to find QN frequencies, and there were many, have 
been devised to compute the lowest QN frequencies (see 
for a review). Secondly, there was no suggestion for 
the asymptotic In 3 value. To know this value serves as 
a strong stimulation. The third reason is obvious: this 
was a very difficult technical problem. 

It is thus satisfying that Motl and Neitzke have proved 
exactly that the limit is In 3. In addition, Van den Brink 
[2^ used a somewhat different approach to rederive the 
value, and at the same time computed the first or- 
der correction to gravitational QN frequencies, confirm- 
ing analytically Nollert's numerical results. However, 
Motl and Neitzke's monodromy method is more flexible 
and powerful: it has been easily adapted to other space- 
times. For example, Castello-Branco and Abdalla 
have generalized the results to the Schwarzschild-de Sit- 
ter geometry, while Schiappa, Natario and Cardoso 
have used it to compute the asymptotic values in the 
-D-dimensional Reissner-Nordstrom geometry. Moreover 
the first order corrections are also easily obtained, as has 
been showed by Musiri and Siopsis .29] for the four di- 
mensional Schwarzschild geometry. Here we shall also 
generalize these first order corrections so that they can 
encompass a general Z?-dimensional Schwarzschild black 
hole. A body of work has been growing on the subject 
of computing highly damped QNMs, both in asymptot- 
ically flat and in de Sitter or anti-de Sitter spacetimes 

As an aside, we note that one of the most intrigu- 
ing features of Motl and Neitzke's technique is that the 
region beyond the event horizon, which never enters in 
the definition of quasinormal modes, plays an extremely 
important role. In particular, the singularity at r = 
is decisive for the computation. At the singularity the 
effective potential for wave propagation blows up. Some- 
how, the equation, or the singularity, knows what we are 
seeking! It is also worth of note the following: for high 
frequencies, or at least frequencies with a large imaginary 
part, the important region is therefore r = where the 
potential blows up, whereas for low frequencies, of inter- 
est for late-time tails [s^, it is the other limit, r — > 00 
which is important. 

The purpose of this paper is two-fold: first we want 
to enlarge and settle the results for the four dimensional 
Schwarzschild black hole. We shall first numerically con- 
firm Nollert's results for the gravitational case, and Berti 
and Kokkotas' 25] results for the scalar case, and see that 
the match between these and the analytical prediction is 
quite good. Then we shall fill-in a gap left behind in these 
two studies: the electromagnetic asymptotic QN frequen- 
cies, for which there are no numerical results, but a lot 
of analytical predictions. This will put the monodromy 
method on a firmer ground. One might say that in the 



3 



present context (of black hole quantization and relation 
with Loop Quantum Gravity) electromagnetic perturba- 
tions are not relevant, only the gravitational ones are. 
This is not true though, because when considering the 
Reissner-Nordstrom case gravitational and electromag- 
netic perturbations are coupled and it is therefore 
of great theoretical interest to study each one separately, 
in the case they do decouple, i.e., the Schwarzschild ge- 
ometry. We find that all the analytical predictions are 
correct: scalar and gravitational QN frequencies asymp- 
tote to 7^ = In 3-|-«(2n-|- l)7r with corrections, where 

-Iff ^ ' V" 

Th = g^fjy is the Hawking temperature. The corrections 
are well described by the existing analytical formulas. On 
the other hand, electromagnetic QN frequencies asymp- 
tote to 7#- = + i2mr with no corrections, as pre- 
dieted by Musiri and Siopsis [29j. Indeed we find that 
the corrections seem to appear only at the level. 

The second aim of this paper is to establish numer- 
ically the validity of the monodromy method by ex- 
tending the numerical calculation to higher dimensional 
Schwarzschild black holes, in particular to the five dimen- 
sional one. We obtain numerically in the five dimensional 
case the asymptotic value of the QN frequencies, which 
turns out to be ln3-|-i(2n-|-l)7r (in units of the black hole 
temperature), for scalar, gravitational tensor and gravi- 
tational vector perturbations. Equally important, is that 
the first order corrections appear as -^^^ for these cases. 
We shall generalize Musiri and Siopsis' method to higher 
dimensions and find that (i) the agreement with the nu- 
merical data is very good (ii) In generic D dimensions 
the corrections appear at the ^(D-3)/{n-2) level. Elec- 
tromagnetic perturbations in five dimensions seem to be 
special: according to the analytical prediction by Motl 
and Neitzke, they should not have any asymptotic QN 
frequencies. Our numerics seem not to go higher enough 
in mode number to prove or disprove this. They indicate 
that the real part approaches zero, as the mode number 
increases, but we have not been able to follow the QN 
frequencies for very high overtone number. We also give 
complete tables for the lowest lying QN frequencies of 
the five dimensional black holes, and confirm previous re- 
sults 0, 35, 36, 37] for the fundamental QN frequencies 
obtained using WKB-like methods. We find no purely 
imaginary QN frequencies for the gravitational QNMs, 
which may be related to the fact that the potentials are 
no longer superpartners. 

II. QN FREQUENCIES OF THE 
FOUR-DIMENSIONAL SCHWARZSCHILD 
BLACK HOLE 

The first computation of QN frequencies of four- 
dimensional Schwarzschild black holes was carried out 
by Chandrasekhar and Detweiler [s^. They were only 
able to compute the lowest lying modes, as it was and 
still is quite hard to find numerically QN frequencies 
with a very large imaginary part. After this pioneer- 



ing work, a number of analytical, semi-analytical and 
numerical techniques were devised with the purpose of 
computing higher order QNMs and to establish the low- 
lying ones using different methods 0, 0| . The most suc- 
cessful numerical method to study QNMs of black holes 
was proposed by Leaver j39| . In a few words his method 
amounts to reducing the problem to a continued frac- 
tion, rather easy to implement numerically. Using this 
method Leaver was able to determine the first 60 modes 
of the Schwarzschild black hole, which at that time was 
about an order of magnitude better than anyone had ever 
gone before. The values of the QN frequencies he deter- 
mined, and we note he only computed gravitational ones, 
suggested that there is an infinite number of QN modes, 
that the real part of the QN frequencies approach a finite 
value, while the imaginary part grows without bound. 
That the real part of the QN frequencies approach a finite 
value was however opened to debate it seemed one 
still had to go higher in mode number in order to ascer- 
tain the true asymptotic behaviour. NoUert showed 
a way out of these problems: he was able to improve 
Leaver's method in a very simple fashion so as to go very 
much higher in mode number: he was able to compute 
more than 2000 modes of the gravitational QNMs. His 
results were clear: the real part of the gravitational QN 
frequencies approaches a constant value ~ 0.0437123/M, 
and this value is independent of I. The imaginary part 
grows without bound as (2?! -f l)/8 where n is any large 
integer. These asymptotic behaviours have both leading 
corrections, which were also determined by NoUert. He 
found that to leading order the asymptotic behaviour of 
the gravitational QN frequencies of a Schwarzschild black 
hole is given by Eq. ^ plus a leading correction, namely, 

wM = 0.0437123+ (2n+ 1) + , (6) 

where the correction term a depends on I. About ten 
years later, Motl proved analytically that this is in- 
deed the correct asymptotic behaviour. To be precise he 
proved that 

= In 3 + i (2n + l)7r -f corrections . (7) 

Th 

Since the Hawking temperature Th is Th = -g^jji 
Nollert's result follows. He also proved that this same re- 
sult also holds for scalar QN frequencies, whereas electro- 
magnetic QN frequencies asymptote to a zero real part. 
Motl and Neitzke [2^ using a completely diS^erent ap- 
proach, have rederived this result whereas Van den Brink 
|26| was able to find the leading asymptotic correction 
term a in JBJ for the_gravitational case. Very recently, 
Musiri and Siopsis 29], leaning on Motl and Neitzke's 
technique, have found the leading correction term for any 
field, scalar, electromagnetic and gravitational. This cor- 
rection term is presented in section Fill Bl where we also 
generalize to higher dimensional black holes. In particu- 
lar, for electromagnetic perturbations they find that the 
first order corrections are zero. 
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FIG. 1: Scalar (s = 0), electromagnetic {s = 1) and 
gravitational {s = 2) QN frequencies of a four-dimensional 
Schwarzschild black hole. Note that the dots indicate the QN 
frequencies there, and the lines connecting the dots only help 
to figure out to which multipole {I, s) they belong to. We show 
the lowest 500 modes for the lowest radiatable multipoles of 
each field, i.e., I = for the scalar field, I — 1 for the electro- 
magnetic and / = 2 for gravitational. However, the asymp- 
totic behaviour is /-independent. It is quite clear from this 
plot that highly damped electromagnetic QN frequencies have 
a completely different behaviour from that of scalar or gravita- 
tional highly damped QN frequencies. In fact, the electromag- 
netic ones asymptote to ^ —> i{2n)n ,n — > oo, as predicted 
by Motl jg^l and Motl and Neitzke ^ . The scalar and grav- 
itational ones asymptote to > In 3 -I- i(2n -I- l)7r , n ^ oo. 



Here we use Nollert's method to compute scalar, elec- 
tromagnetic and gravitational QN frequencies of the four- 
dimensional Schwarzschild black hole. The details on the 
implementation of this method in four dimensions are 
well known and we shall not dwell on it here any more. 
We refer the reader to the original references |2lL l39| . 
Our results are summarized in Figures Ill2l and Tables 
m In Fig. n we show the behaviour of the 500 low- 
est QN frequencies for the lowest radiatable multipole of 
each field (we actually have computed the first 5000 QN 
frequencies, but for a better visualization we do not show 
them all): ^ = for scalar, Z = 1 for electromagnetic and 
I = 2 gravitational. Note that in Fig. as well as in 
all other plots in this work, the dots indicate the QN 
frequencies there, and the lines connecting the dots only 
help to figure out to which multipole (I, s) they belong 
to. In what concerns the asymptotic behaviour, our nu- 
merics show it is not dependent on I, and that formula 
© holds. The results for the scalar and gravitational QN 
frequencies are not new: the gravitational ones have been 
obtained by NoUert , as we previously remarked, and 
the scalar ones have been recently arrived at by Berti and 
Kokkotas 25,]. Here we have verified both. Both scalar 
and gravitational QN frequencies have a real part asymp- 



TABLE L The correction coefficients for the four-dimensional 
Schwarzschild black hole, both numerical, here labeled as 
"Corr4" and analytical, labeled as "Corr4^". These results 
refer to scalar perturbations. The analytical results are ex- 
tracted from |2g|| (see also formula 1291 below *). Notice the 
very good agreement between the numerically extracted re- 
sults and the analytical prediction. These values were obtained 
using the first 5000 modes. 



s = 


I 


Corr^ 


Corr^ 





1.20-1.20i 


1.2190-1.21901 


1 


8.52- 8.52i 


8.5332-8.53321 


2 


23.15-23.15i 


23.1614-23.16141 


3 


44.99-44.99i 


45.1039-45.10391 


4 


74.34-74.34i 


74.3604-74.36051 





toting to In 3 X Th- Electromagnetic QN frequencies, on 
the other hand, have a real part asymptoting to 0. This 
is clearly seen in Fig. ^ For gravitational perturbations 
Nollert found, and we confirm, that the correction term 
a in © is 

a = 0.4850, I = 2 
a = 1.067, 1 ^3 (8) 

We have also computed the correction terms for the scalar 
case, and we have obtained results in agreement with 
Berti and Kokkotas' ones [2^ . Our results are summa- 
rized in Tables IIIIII where we also show the analytical 
prediction (see section UlIBII . To fix the conventions 
adopted throughout the rest of the paper, we write 

w , , N Corr 

— =ln3 + 2(2n+l)7r+^^. (9) 

Thus, the correction term a in Q is related to Corr by 
a = Corr x Th . 

As for the electromagnetic correction terms, we found 
it was not very easy to determine them. In fact, it is 
hard, even using Nollert's method, to go very high in 
mode number for these perturbations (we have made it 
to n = 5000). The reason is tied to the fact that these 
frequencies asymptote to zero, and it is hard to deter- 
mine QN frequencies with a vanishingly small real part. 
Nevertheless, there are some features one can be sure of. 
We have fitted the electromagnetic data to the following 
form — i(2n)n + , and found that this gave very 
poor results. Thus one saw numerically that the first or- 
der correction term is absent, as predicted in '2^ (see also 
section Fill Bl where we rederive these corrections, gener- 
alizing them to arbitrary dimension). Our data seems to 
indicate that the leading correction term is of the form 
^j^. This is more clearly seen in Fig. [21 where we plot 
the first electromagnetic QN frequencies on a In plot. 
For frequencies with a large imaginary part, the slope is 
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TABLE II: The correction coefficients for the four- 
dimensional Schwarzschild black hole, both numerical, here 
labeled as "Corr4 " and analytical, labeled as "Corr4^". These 
results refer to gravitation al p erturbations. The analytical 
results are extracted from |2fll | (see also formula l|29l below). 
Notice the very good agreement between the numerically ex- 
tracted results and the analytical prediction. These values 
were obtained using the first 5000 modes. 
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Corr^: 


Corr^ 
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6.08-6.081 


6.0951-6.09511 
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13.39-13.39i 


13.4092-13.40921 
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23. 14-23.141 


23.1614-23.16141 


5 


35.33-35.33i 


35.3517-35.35171 


6 


48.90-48.90i 


49.9799-49.97991 
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FIG. 2: Electromagnetic QN frequencies of a four-dimensional 
Schwarzschild black hole. 



about —2/3 as it should if the asymptotic behaviour is 
of the form 



LO b 

Iff n-^/^ 



(10) 



Our numerics indicate that b has an Z-dependence going 
like 



b ^ Cil{l + 1) + C2 [l{l + l)f + C3 [l{l + l)f , (11) 



where ci, C2, andcs are constants. This is also the ex- 
pected behaviour, since a -^^j2 dependence means we have 
to go to third order perturbation theory, where we get 
corrections of the form Hll|) as one can easily convince 
oneself. 



III. QN FREQUENCIES OF D-DIMENSIONAL 
SCHWARZSCHILD BLACK HOLES 



A. Equations and conventions 

The metric of the D-dimcnsional Schwarzschild black 
hole in {t, r, 6i, 6*2, .., 6*0-2) coordinates is [ill 11^ 



ds^ = -fdf + r^dr^ + r^dn^ 



1-2 , (12) 



with 



f-^-Zn^- (13) 



The mass of the black hole is given by M ~ 



16ng ' 

where ^Id~2 — r [ {D-i)/2] area of a unit (_D — 2) 

sphere, (if2|,_2 the line element on the unit sphere 
S^~'^, and Q is Newton's constant in £)-dimensions. In 
the following it will prove more useful to rescale vari- 
ables so that the form of the metric is l|12|l but with 
/ = 1 — j.n-3 ) i-e., we shall choose m — 1, following the 
current fashion. We will only consider the linearized ap- 
proximation, which means that we are considering wave 
fields outside this geometry that are so weak they do not 
alter this background. The evolution equation for a mass- 
less scalar field follows directly from the Klein-Gordon 
equation (see _43J for details). The gravitational evolu- 
tion equations have recently been derived by Kodama 
and Ishibashi . There are three kinds of gravitational 
perturbations, according to Kodama and Ishibashi's ter- 
minology: the scalar gravitational, the vector gravita- 
tional and the tensor gravitational perturbations. The 
first two already have their counterparts in Z? = 4, which 
were first derived by Regge and Wheeler and by Zer- 
illi [4^ . The tensor type is a new kind appearing in higher 
dimensions. However, it obeys exactly the same equation 
as massless scalar fields, as can easily be seen. Due to 
the complex form of the gravitational scalar potential, 
we shall not deal with it. Instead, we shall only consider 
the tensor and vector type of gravitational perturbations. 
In any case, if the analytic results are correct, then the 
gravitational scalar QN frequencies should have the same 
asymptotic form as the gravitational vector and tensor 
QN frequencies. The evolution equation for the electro- 
magnetic field in the higher dimensional Schwarzschild 
geometry was arrived at for the first time by Crispino, 
Higuchi and Matsas It has recently been rederived 
by Kodama and Ishibashi, in a wider context of charged 
black hole perturbations. We shall follow Kodama and 
Ishibashi's terminology. According to them, there are 
two kinds of electromagnetic perturbations: the vector 
and scalar type. If one makes the charge of the black hole 
Q = in Kodama and Ishibashi's equations one recovers 
the equations by Crispino, Higuchi and Matsas, although 
this seems to have been overlooked in the literature. The 
evolution equation for all kinds of fields (scalar, gravita- 
tional and electromagnetic) can be reduced to the second 
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order differential equation 



(14) 



where r is a function of the tortoise coordinate , defined 
through — f{r), and the potential can be written in 
compact form as 



V = f{r) 



l{l + D-3) iD-2){D-A) 



4^2 



+ 



21 



47-D-l 



(15) 



The constant j depends on what kind of field one is study- 
ing: 



0, 



2 

D-2 



D-2 ' 



scalar and gravitational tensor 
perturbations, 
gravitational vector 
perturbations, 
electromagnetic vector 
perturbations, 
electromagnetic scalar 
perturbations . 



(16) 



Notice that in four dimensions , j r educes to the usual val- 
ues given by Motl and Neitzke |23| . According to our con- 
ventions, the Hawking temperature of a Z)-dimensional 
Schwarzschild black hole is 



H 



D-3 
47r 



(17) 



Our purpose here is to investigate the QN frequencies 
of higher dimensional black holes. Of course one cannot 
study every D. We shall therefore focus on one particu- 
lar dimension, five, and make a complete analysis of its 
QN frequencies. The results should be representative. 
In particular, one feature that distinguishes an arbitrary 
D from the four dimensional case is that the correction 
terms come with a different power, i.e., whereas in four 
dimensions the correction is of the form we shall find 

that in generic D it is of the form ^(D-3)/(D-2-j ■ Thus, if 
one verifies this for D = 5 for example, one can ascertain 
it will hold for arbitrary D. We shall now briefly sketch 
our numerical procedure for finding QN frequencies of 
a five dimensional black hole. This is just Leaver's and 
Nollert's technique with minor modifications. 

The QNMs of the higher-dimensional Schwarzschild 
black hole are characterized by the boundary conditions 
of incoming waves at the black hole horizon and outgoing 
waves at spatial infinity, written as 



(18) 



where the time dependence of perturbations has been 
assumed to be e*"^*. 



B. Perturbative calculation of QNMs of 
D-dimensional Schwarzschild black holes 

In this section we shall outline the procedure for com- 
puting the first order corrections to the asymptotic value 
of the QN frequencies of a D-dimensional Schwarzschild 
black hole. This will be a generalization of Musiri and 
Siopsis' method 0|, so we adopt all of their notation, 
and we refer the reader to their paper for further details. 
The computations are however rather tedious and the fi- 
nal expressions are too cumbersome, so we shall refrain 
from giving explicit expressions for the final result. 

We start with the expansion of the potential V near 
the singularity r — 0. One can easily show that near this 
point the potential (|15|l may be approximated as 



V • 



1-/ 



A 



(-Zw)(^-3)/(-D-2) 



(19) 



where we adopted the conventions in [22j and therefore 
z = LOTtf. The constant A is given by 



A ^ 



+ D-3) + {D-2){D-A)-{D- 2)^(1 - f) 

[D -2YD-l)/{D-2) 

2(l_j2)(^_2)(2C-5)/(^-2) 



2D 



(20) 



In D = 4 this reduces to the usual expression [23, |43| 
for the potential near r = 0. Expression (|19|l is a formal 
expansion in ^(d-3]/{d-2) i so we may anticipate that in- 
deed the first order corrections will appear in the form 
„(-D-3)/(-D-2) ■ So now we may proceed in a direct manner: 
we expand the wavefunction to first order in ^(d-3]/{d-2) 
as 



* = * 



1 



.<5(i) 



(21) 



^(D-3)/(D-2) 

and find that the first-order correction obeys the equation 
dz^ Az^ 



with 



5V = - 



A 



^(D-3)/(I5-2)(_^)(I5-l)/(D-2) 



(23) 



All of Musiri and Siopsis' expressions follow directly to 
the Z?-dimensional case, if one makes the replacement 
^/—LorQ ijj{D-'i)/(D-2) jj-^ their expressions. For ex- 
ample, the generalsolution of H22|l is 



where 



(0) 



.(D-3)/(Z5-2) 



sin j'7r/2 



(24) 
(25) 



7 



and the wavefunctions are 



(0) 



(26) 



The only minor modification is their formula (30). In 
U-dimensions, it follows that 



G±(z), 



(27) 



where k = 2(0-2) ' ^"^^ ^'^^ even analytic functions 
of z. So we have all the ingredients to construct the first 
order corrections in the D-dimensional Schwarzschild ge- 
ometry. Unfortunately the final expressions turn out to 
be quite cumbersome, and we have not managed to sim- 
plify them. We have worked with the symbolic manip- 
ulator Mathematica. It is possible to obtain simple ex- 
pression for any particular D, but apparently not for a 
generic D. For generic D we write 

^=ln(l + 2cos^j)+z(2n+l)7r + --p3^^^-^ (28) 

So the leading term is ^ = In 3 + i{2n + l)7r, for scalar, 
gravitational tensor ancf gravitational tensor perturba- 
tions (and the same holds for gravitational scalar) . How- 
ever, for electromagnetic perturbations (vector or scalar) 
in D ~ 5, 1 + 2 cos irj — 1 + 2 cos ^ = 0. So there seem 
to be no QN frequencies as argued by Motl and Neitzke 
[2^ . We obtained for D — A the following result for the 
correction coefficient in (I28f) . 



Corr4 



(1 - i)7r3/2(_i + f - 31(1 + i))££i2£r[l/4] 



3(1 + 2 cos j7r)r[3/4]r[3/4 - j72]r[3/4 + j/2] ' 

(29) 

which reduces to Musiri and Siopsis' expression. For D = 
5 it is possible also to find a simple expression for the 
corrections: 



Corr5 — 



^(9J•'-(24 + 20^(^ + 2)))7^3/2 
120 X 3i/3r[2/3]r[2/3 - i/2]r[2/3 + j/2] 



;;rfl/6l "'"»V3+j/2)7r) sin((2/3+j/2)7r) 

The limits j ^ 0, 2 are well defined and yield 

^ ,=0 (l-^^/3)(6 + 10^ + 5/2)7r3/2^[l/6] 
Corrs'' = 



Corrs^^^ = 



45 X 3i/3r[2/3]3 

(-1 + iV3){-3 + 10/ + 5P)7r3/2r[l/6] 
45 X 3i/3r[-l/3]r[2/3]r[5/3] 



(31) 



_ (32) 
We list in Tables IVIIVIII some values of this five dimen- 
sional correction for some values of the parameter j and I 
and compare them with the results extracted numerically. 
The agreement is quite good. The code for extracting the 
generic £)-dimensional corrections is available from the 
authors upon request. It should not come as a surprise 
that for j — 2/3 the correction term blows up: indeed 
already the zeroth order term is not well defined. 



C. Numerical procedure and results 

1. Numerical procedure 

In order to numerically obtain the QN frequencies, 
in the present investigation, we make use of NoUert's 
method |^, since asymptotic behaviors of QNMs in the 
limit of large imaginary frequencies are prime concern in 
the present study. 

In our numerical study, we only consider the QNMs 
of the Schwarzschild black hole in the five dimensional 
spacetime, namely the D = 5 case. As we have remarked, 
this should be representative. The tortoise coordinate is 
then reduced to 



1 



where x ~ 



+ — ln(a; - xi) 
2x1 



and xi 



1 

2^ 



ln(x-|-a;i), (33) 



Here, rh stands for the 



horizon radius of the black hole. The perturbation func- 
tion ^' may be expanded around the horizon as 



* = e 



00 

\x-x,rix+x^ryaj 



k=0 



-Xi 



(34) 



where p — iuj/2xi and ao is taken to be ao — 1. The 
expansion coefficients in equation 134|l are determined 
via the four-term recurrence relation (it's just a matter of 
substituting expression in the wave equation l(T^ '). 
given by 

aofli + /3oao = , 

aia2 + /3iai + 7iao = , (35) 
ako-k+i + Pkak + Ikdk-i + Skak-2 = , fc = 2, 3, • • • , 



where 



Pk 

Ik 
6k 



2(2p+fc + l)(fc+l), 

-5(2p + k){2p + k + l)~l{l + 2) 



4{2p+k 



l)(2p + fc+l) + -(l-/), 
-{2p +k- 2){2p + k + l)-^{l-f). 



It is seen that since the asymptotic form of the pertur- 
bations as r* ^ c» is written in terms of the variable x 
as 



[x - xiy{x + xiY , 



(36) 



the expanded perturbation function 5" defined by equa- 
tion (|34(l automatically satisfy the QNM boundary con- 
ditions (|18|l if the power series converges for < a; < Xi. 
Making use of a Gaussian elimination l49| , we can reduce 
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the four-term recurrence relation to the three-term one, 
given by 

a'oai + PqUq = , 

a'^.ak+l + P'kak + TfeOfc-i = , fc = 1, 2, • • • , (37) 

where aj,, /3^, and 7^ are given in terms of a^, (3^, 7^ and 
h by 



A = P'k = f^k, 7fe = Ik, for fc = 0, 1, (38) 



and 



O^k = Oik, 

P'k = Pk- a'k^iSk/l'k^i , 

Ik = Ik- P'k-i^k/lk-i^ forfc>2. 



(39) 



Now that we have the three-term recurrence relation for 
determining the expansion coefficients ak, the conver- 
gence condition for the expansion 134|) , namely the QNM 
conditions, can be written in terms of the continued frac- 
tion as |39.l50l| 



O^'oll «l72 0^273 

P'l- P's" ■ 



"o7i 



0,(40) 



TABLE III: The first lowest QN frequencies for scalar and 
gravitational tensor (j — 0) perturbations of the five dimen- 
sional Schwarzschild black hole. The frequencies are normal- 
ized in units of black hole temperature, so the Table really 
shows 7^, where Th is the Hawking temperature of the black 
hole. 



J = 


n 


/ = 0: 


I = 1: 


I = 2 





3.3539+2.40891 


6.3837+2.2764i 


9.4914+2. 2462i 


1 


2.3367-^8.31011 


5.3809+7.2734i 


8.7506+6.9404i 


2 


1.8868+14.7861 


4.1683+13. 252i 


7.5009+12. 225i 


3 


1.6927+21.2191 


3.4011+19.708i 


6. 2479+18. 214i 


4 


1.5839+27.6071 


2. 9544+26. 215i 


5.3149+24.597i 



fl'2- 



2. Numerical results 

Using the numerical technique described in section 
IIII C II we have extracted the 5000 lowest lying QN fre- 
quencies for the five dimensional Schwarzschild black 
hole, in the case of scalar, gravitational tensor and grav- 
itational vector perturbations. For electromagnetic vec- 
tor perturbations the situation is different: the real part 
rapidly approaches zero, and we have not been able to 
compute more than the first 30 modes. 



where the first equality is a notational definition com- 
monly used in the literature for infinite continued frac- 
tions. Here we shall adopt such a convention. In order 
to use Nollert's method, with which relatively high-order 
QNM with large imaginary frequencies can be obtained, 
we have to know the asymptotic behaviors of ak+i/ak in 
the limit of fc — > cxd. According to a similar consideration 
as that by Leaver j42j|, it is found that 



Clfc-l-l 

ak 



(41) 



where the sign for the second term in the right-hand side 
is chosen so as to be 



Re(±2^) < . 



(42) 



In actual numerical computations, it is convenient to 
solve the fc-th inversion of the continue fraction equation 
gnj), given by 



«fc-i7fc «fc_27fc- 

l^'k-l~ l^k-2~ 

l^'k+l^ P'k+2^ 



ao7i 



(43) 



The asymptotic form (|41ll plays an important role in 
Nollert's method when the infinite continued fraction in 
the right-hand side of equation 14311 is evaluated 



(i) Low-lying modes 

Low lying modes are important since they govern the 
intermediate-time evolution of any black hole perturba- 
tion. As such they may play a role in TeV-scale gravity 
scenarios Wa. and higher dimensional black hole forma- 
tion [33.l52l |. For example, it is known 0111,13 that if 
one forms black holes through the high energy collision 
of particles, then the fundamental quasinormal frequen- 
cies serve effectively as a cuttof in the energy spectra of 
the gravitational energy radiated away. In Tables IIIIIVI 
we list the five lowest lying QN frequencies for some val- 
ues of the multipole index The fundamental modes 
are in excellent agreement with the ones presented by 
Konoplya [H, [s^ using a high-order WKB approach. 
Notice he uses a different convention so one has to be 
careful when comparing the results. For example, in our 
units Konoplya obtains for scalar and tensorial pertur- 
bations (j = 0) with / = 2 a fundamental QN frequency 
ujq/Th = 9.49089 + 2.24721i, whereas we get, from Ta- 
ble Uni the number 9.4914 + 2.2462i so the WKB ap- 
proach does indeed yield good results, at least for low- 
lying modes, since it is known it fails for high-order ones. 

In the four dimensional case, and for gravitational vec- 
tor perturbations, there is for is each Z, a purely imag- 
inary QN frequency, which had been coined an "alge- 
braically special frequency" by Chandrasekhar |33j. For 
further properties of this special frequencies we refer the 
reader to [55j| . The existence of these purely imagi- 
nary frequencies translates, in the four dimensional case. 
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TABLE IV: The first lowest QN frequencies for gravita- 
tional vector {j — 2) perturbations of the five dimensional 
Schwarzschild black hole. The frequencies are normalized in 
units of black hole temperature, so the Table really shows 7^ , 
where Th is the Hawking temperature of the black hole. 



500 



J = 2 


n 


I = 2: 


I = 3: 


/ = 4 





7.1251+2.05791 


10.8408+2.09761 


14.3287+2.13641 


1 


5.9528+6.42171 


8.8819+1.09291 


12.8506+1.05741 


2 


3.4113+12.09291 


7.2160+16.39791 


11.5069+16.02741 


3 


2.7375+19.60941 


5.5278+22.33211 


9.99894+21.50091 


4 


2.5106+26.26251 


3.9426+28.65691 


8.56386+27.43391 



TABLE V; The first lowest QN frequencies for electromag- 
netic vector (j = 2/3) perturbations of the five dimensional 
Schwarzschild black hole. The frequencies are normalized in 
units of black hole temperature, so the Table really shows 7^ , 
where Th is the Hawking temperature of the black hole. 



j = 2/3 


n 


I = 1: 


I = 2: 


/ = 3 





5.9862+2.20381 


9.2271+2.21441 


12.4184+2.21771 


1 


4.9360+7.06761 


8.4728+6.84591 


11.8412+6.76601 


2 


3.6588+12.95811 


7.1966+12.08041 


10.7694+11.66261 


3 


2.8362+19.34221 


5.9190+18.03431 


9.4316+17.10191 


4 


2.3322+25.78611 


4.9682+24.38321 


8.1521+23.07691 



a relation between the two gravitational wavefunctions 
(i.e., between the Regge- Wheeler and the Zerilli wave- 
function, or between the gravitational vector and gravi- 
tational scalar wavefunctions, respectively). It is possi- 
ble to show for example that the associated potentials 
are related through supersymmetry. Among other con- 
sequences, this relation allows one to prove that the QN 
frequencies of both potentials are exactly the same [23| • 
We have not spotted any purely imaginary QN fre- 
quency for this five-dimensional black hole. In fact the 
QN frequency with the lowest real part is a gravitational 
vector QN frequency with I — 4 and overtone number 
n = 13, w = 8.7560 x IQ-^ + 6.9934i. This may indicate 
that in five dimensions, there is no relation between the 
wavefunctions, or put another way, that the potentials 
are no longer superpartner potentials. This was in fact 
already observed by Kodama and Ishibashi 0| for any 
dimension greater than four. It translates also in Tables 
IIIIIIVI which yield different values for the QN frequen- 
cies. Although we have not worked out the gravitational 
scalar QN frequencies, a WKB approach can be used for 
the low-lying QN frequencies, and also yields different 
values pja . l3a | . 

(ii) Highly damped modes 
Our numerical results for the highly damped modes. 




2 4 6 



FIG. 3: The QN frequencies of scalar (j — 0), gravitational 
vector (j = 2) and electromagnetic vector (j = 2/3) pertur- 
bations of a five-dimensional Schwarzschild black hole. We 
only show the lowest radiatable multipoles, but the asymtp- 
toic behaviour is /-independent. 



i.e., QN frequencies with a very large imaginary part, 
are summarized in Figures I3I5I and in Tables IVIIVIII 

In Fig. 131 we show our results for scalar, gravitational 
tensor (j = 0), gravitational vector {j = 2) and electro- 
magnetic vector (j = 2/3) QN frequencies. For j = 0,2 
the real part approaches In 3 in units of the black hole 
temperature, a result consistent with the analytical re- 
sult in 23J (see also section IlII B|l . Electromagnetic vec- 
tor {j = 2/3) QN frequencies behave differently: the real 
part rapidly approaches zero, and this makes it very dif- 
ficult to compute the higher modes. In fact, the same 
kind of problem appears in the Kerr geometry |3C| . and 
this a major obstacle to a definitive numerical charac- 
terization of the highly damped QNMs in this geometry. 
We have only been able to compute the first 30 modes 
with accuracy. The prediction of ,23|] for this case (see 
section lTlI B|l is that there are no asymptotic QN frequen- 
cies. It is left as an open question whether this is true or 
not. Our results indicate that the real part rapidly ap- 
proaches zero, but we cannot say whether the modes die 
there or not, or even if they perform some kind of oscil- 
lation. The imaginary part behaves in the usal manner, 
growing linearly with mode number. 

Let us now take a more detailed look at our numeri- 
cal data for scalar and gravitational QN frequencies. In 
Fig. ^ we show in a In plot our results for scalar and 
gravitational tensor (j = 0) QN frequencies. 

Our numerical results are very clear: asymptotically 
the QN frequencies for scalar and gravitational tensor 
perturbations (j = 0) behave as 



UJ 



ln3 + i (2n + l)n + 



Corrs 



(44) 
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0.01 0.1 1 10 100 0.001 0.01 0.1 1 10 100 

Re(a;/TH)-ln3 Re{u/T^^)-ln3 



FIG. 4: The QN frequencies of scalar and gravitational ten- 
sor perturbations {j = 0) of a five- dimensional Schwarzschild 
black hole. Notice that for frequencies with a large imaginary 
part the slope is /-independent and equals 3/2. 



TABLE VI: The correction coefficients for the five- 
dimensional Schwarzschild black hole, both numerical, here 
labeled as "Corr^" and analytical, labeled as "Corr^". These 
results refer to scalar or gravitational tensor perturbations 
{j — 0). The analytical results are extracted from the analyt- 
ical formula 1311 . Notice the very good agreement, to within 
0.5% or less, between the numerically extracted results and 
the analytical prediction. 



j = 



I 


Corr^: 




Corrg" 





1.155 - 


1.993i 


1.15404-1.998861 


1 


4.058 - 


6.956i 


4.03915-6.996011 


2 


8.894 - 


15.26i 


8.84765-15.32461 


3 


15.64 - 


26.79i 


15.5796-26.98471 


4 


24.33 - 


41.67i 


24.2349-41.97611 



So the leading terms In 3 -I- (2n -I- l)i are indeed the ones 
predicted in ,23j. Interestingly, the first corrections do 
not appear as but as :^^t73- This was shown to be the 
expected analytical result in section Fill Bl where we gen- 
eralized Musiri and Siopsis' (29j results to higher dimen- 
sions. In table lVll we show the coefficient Corr5 extracted 
numerically along with the predicted coefficient (see ex- 
pression pifl ). The table is very clear: the numerical 
values match the analytical ones. Another confirmation 
that the corrections appear as ^^^73- is provided by Fig. 01 
In this figure the QN frequencies are plotted in a In plot 
for an easier interpretation. One sees that for large imag- 
inary parts of the QN frequencies, the slope of the plot 
is approximately —3/2, as it should be if the corrections 
are of the order -^r^ . 

In Fig. [SI we show in a In plot our results for grav- 



FIG. 5: The QN frequencies of gravitational vector perturba- 
tions {j = 2) of a five- dimensional Schwarzschild black hole. 



TABLE VII: The correction coefficients for the five- 
dimensional Schwarzschild black hole, both numerical, here 
labeled as "Corr^" and analytical, labeled as "Corr^", for 
gravitational vector perturbations {j = 2). The analytical re- 
sults are extracted from the analytical formula 1321 . Notice 
the very good agreement, to within 0.5% or less, between the 
numerically extracted results and the analytical prediction. 



1=1 



I 


Corr^: 




Corr^: 


2 


3.568 - 


6.1391 


3.5583 -6.163161 


3 


6.947 - 


11.941 


6.92425-11.99321 


4 


11.29 - 


19.38 i 


11.2519-19.48891 



itational vector {j — 2) QN frequencies. Again, vector 
QN frequencies have the asymptotic behaviour given by 
expression ()44|l . with a different correction term Corrs. 
Again, the corrections show themselves in the In 
plot of Fig. ^ for very large imaginary parts, the slope 
is —3/2, as it should. The numerically extracted coeffi- 
cient Corrs for j = 2 perturbations is hsted in Table IVIll 
along with the analytically predicted value (see section 
llllljll . the agreement is very good. 



IV. DISCUSSION OF RESULTS AND FUTURE 
DIRECTIONS 

We have made an extensive survey of the QNMs of 
the four and five dimensional Schwarzschild black hole. 
The investigation presented here makes a more com- 
plete characterization of the highly damped QNMs in 
the Schwarzschild geometry. In the four-dimensional 
case, we confirmed previous numerical results regarding 
the scalar and gravitational asymptotic QN frequencies. 
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We found that both the leading behaviour and the first 
order corrections for the scalar and gravitational per- 
turbations agree extremely well with existing analyti- 
cal formulas. We have presented new numerical results 
concerning electromagnetic QN frequencies of the four- 
dimensional Schwarzschild geometry. In particular, this 
is the first work dealing with highly damped electromag- 
netic QNMs. Again, we find that the leading behaviour 
and the first order corrections agree with the analytical 
calculations. The first order corrections appear at the 
^jj72- level. In the five-dimensional case this represents 
the first study on highly damped QNMs. We have seen 
numerically that the asymptotic behaviour is very well 
described by ^ = ln3 + TTi{2n +1) for scalar and gravi- 
tational perturbations, and agrees with the predicted for- 
mula. Moreover the first order corrections appear at the 
Isvel which is also very well described by the ana- 

lytical calculations, providing one more consistency check 
on the theoretical framework. In generic D dimensions 
the corrections appear at the ^(o-3)/(d-2) level. We have 
not been able to prove or disprove the analytic result 
for electromagnetic QN frequencies in five dimensions 
(j = 2/3). Other conclusions that can be taken from our 
work are: the monodromy method by Motl and Neitzke 
is correct. It is important to test it numerically, as we 
did, since there are some ambiguous assumptions in that 
method. Moreover, their method is highly flexible, since 
it allows an easy computation of the correction terms, 
as shown by Musiri and Siopsis, and generalized here 
for the higher dimensional Schwarzschild black hole. We 
have basically showed numerically that the monodromy 
method by Motl and Neitzke [23|, and its extension by 
Musiri and Siopsis to include for correction terms in 
overtone number, are excellent techniques to investigate 
the highly damped QNMs of black holes. 

Highly damped QNMs of black holes are the bed rocks 
that the recent conjectures [13 . IT^ 0^ a-re built on, re- 
lating these to black hole area quantization. However, to 
put on solid ground the conjecture that QNMs can ac- 
tually be of any use to black hole area quantization, one 
has to do much better then one has been able to do up 
to now. In particular, it is crucial to have a deeper un- 
derstanding of the Kerr and Reissner-Nordstrom QNMs. 



These are, next to the Schwarzschild, the simpler asymp- 
totically flat geometries. If these conjectures hold, they 
must to do also for these spacetimes. Despite the fact 
that there are convincing numerical results for the Kerr 
geometry ,30j , there seems to be for the moment no seri- 
ous analytical investigation of the highly damped QNMs 
for this spacetime. Moreover, even though we are in pos- 
session of an analytical formula for the hi ghly damped 
QNMs of the Reissner-Nordstrom geometry , and this 
has been numerically tested already , we have no idea 
what it means! In particular, can one use it to quantize 
the black hole area? In this case the electromagnetic per- 
turbations are coupled to the gravitational ones. How do 
we proceed with the analysis, that was so simple in the 
Schwarzschild geometry? Is the assumption of equally 
spaced eigenvalues a correct one, in this case? Can these 
ideas on black hole area quantization be translated to 
non-asymptotically flat spacetimes, as the de Sitter or 
anti-de Sitter spacetime? These are the fundamental is- 
sues that remain to be solved in this field. Should any 
of them be solved satisfactorily, then these conjectures 
would gain a whole new meaning. As they stand, it may 
just be a numerical coincidence that the real part of the 
QN frequency goes to In 3. 
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